Penetration of Andreev bound states into the ferromagnet in a 
SrRuO3/(110)YBa2Cu3O7^ bilayer: a scanning tunneling spectroscopy study 
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Scanning tunneling spectroscopy of thin epitaxial SrRuOz/{)~^Q)y Ba2Cu307-s ferromag- 
net /superconductor bilayers, reveal a clear penetration of the Andreev bound states into the fer- 
romagnetic layer. The penetration is manifested in the density of states of the ferromagnet as a 
split zero bias conductance peak with an imbalance between peak heights. Our data indicate that 
the splitting occurs at the superconductor side as a consequence of induced magnetization, con- 
firming recent theoretical predictions. The imbalance is attributed to the spin polarization in the 
ferromagnet. 

PACS numbers: 74.45.-l-c, 74.50.-fr, 74.78.Fk, 74.81.-g 



I. INTRODUCTION 



The study of superconductor (S)/ ferromagnet (F) 
proximity systems allows a direct investigation of the in- 
terplay between the two competing orders of supercon- 
ductivity and ferromagnetism. In a highly transparent 
S/normal- metal (N) junction, superconducting correla- 
tions are induced in N while they are weakened in the 
S side. The mechanism underlying the proximity effect 
(PE) is the Andreev reflection (AR): a hole- like quasi- 
particle impinging on the interface from the N side is 
retro-reflected as an electron-like quasiparticle with op- 
posite momentum and spin while destroying a Cooper 
pair in the S. Consequently, in S/F junctions, the prop- 
erties of the PE are significantly modified due to the 
spin polarization and the presence of an exchange field 
in F. Theoretical works predict a rapid oscillatory decay 
of the induced superconducting order parameter inside 
Fiii^i^ These predictions, that were confirmed by vari- 
ous experiments^'^'^, stem from the singlet pairing in the 
S, and hold for both s-wave and d-wave (along antinodal 
directions) superconductorsi^*^ However, the understand- 
ing of how the anisotropy of the d-wave symmetry within 
the a-b plane manifests itself in the PE in the presence of 
an exchange field, is still rudimentary. In particular, it is 
still unclear what happens to the Andreev bound states 
(ABS) that reside at the nodal surfaces of d-wave S and 
how, in turn, they affect the density of states (DOS) of 
F. 

Splitting the current carrying ABS requires the re- 
moval of either their directional or spin degeneracy. The 
former may result from the admixture of a subdominant 
order parameter— while the latter may be due to a small 
magnetization at the surface. A possible origin of such 
a surface magnetization is the background antiferromag- 
netic correlations in the underdoped regime of the high- 
Tc superconductors.— Alternatively, a magnetization in- 
side the S may also result from a proximity to a mag- 
netic layer. Recent theoretical works that considered the 
inverse PE in S/F bilayers predict an induced magneti- 
zation in the S side over a length-scale of the supercon- 



ducting coherence length, ^g .^-'^d^d^d^ xhe direction of 
the induced magnetic moment near the interface may be 
parallel or antiparallel to the magnetization in F, depend- 
ing on the interface transparency, the thickness of F and 
the strength of the exchange field. Experimentally, there 
is little evidence for such an induced magnetization] -'^^d^ 
The induced magnetization is expected to remove the 
spin degeneracy of the ABS at the S/F interface and 
shift them to finite energies. However, this effect was not 
observed so far in the DOS of F and is the focus of this 
paper. 

We have previously shown in N/nodal-oriented rf-wave 
S bilayers that as long as the interface is highly transpar- 
ent, the N layer is ballistic and phase coherence is main- 
tained, ABS can penetrate into the N layer and give rise 
to a zero-bias conductance peak (ZBCP) in its DOS.— 
The question that arises now is what happens when an 
exchange field is introduced in the normal layer. The 
possible penetration of the (spin degenerate) ABS into a 
spin polarized F layer was not observed experimentally 
and remains unclear theoretically. 

In this work we employed scanning tunneling spec- 
troscopy on thin epitaxial SrRuO^/ {11Q)Y Ba2Cu'iO'j-s 
(SRO/(110)YBCO) S/F bilayers below fufl coverage of 
the YBCO by the SRO layer. Our spectroscopy mea- 
surements clearly reveal an intricate penetration of the 
ABS into the F layer to distances of at least 9 nm. The 
penetration is refiected in the DOS of the F layer as a 
split ZBCP with a pronounced imbalance between split 
peak heights. The split is attributed to the removal of 
the spin degeneracy of the ABS due to induced magneti- 
zation in the S. The imbalance between the peak heights 
is attributed to the spin polarization in the F layer. 



II. EXPERIMENT 

SRO is an itinerant ferromagnet {Tcurie — in 
thin films) with lattice parameters that are similar to 
those of YBCOiS and therefore they can form highly 
transparent interfaces, essential for the existence of AR 
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FIG. 1: (a) SEM image (0.8X0.6 /im^) of a 3 nm thick SRO 
layer on a (IIO)YBCO film demonstrating the SRO island 
(bright areas) topography. Inset: a zoom on one SRO crys- 
tallite, (b) Magnetization measurement (at 1000 Oe after 
zero-field cooling) depicting the ferromagnetic transition on- 
set at 150 K and the S transition at 88 K. 



and the PE. A total of 16 bilayers of SRO (2-17 nm nom- 
inal thickness) on (110) YBCO were prepared by laser 
ablation deposition on (110) SrTiO^ substrates. In order 
to achieve YBCO layers with the nodal (110) orientation, 
first a 10 nm thick template layer of YBCO was deposited 
at 660°C, and then another YBCO layer of 50 nm thick- 
ness was deposited at 750°C. The (110) orientation was 
verified by x-ray diffraction. Finally, the SRO layer was 
deposited at 785° C. Annealing for obtaining optimally 
doped YBCO was done in-situ under 50 Torr oxygen gas 
pressure and a dwell of one hour at 450° C (the SRO layer 
is insensitive to the oxygen annealing process). The bi- 
layers showed S transition temperatures at around 88 K 
with a transition width of about 3 K, implying nearly 
optimally doped films. The tunneling spectra (dl/dV vs. 
V characteristics, proportional to the local quasi-particle 
DOS) were obtained at 4.2 K, much lower than both the 
superconducting and ferromagnetic transitions. 



III. RESULTS AND DISCUSSION 

The bare (IIO)YBCO films feature - 40 x 100 nm^ 
elongated crystallite with uniform directionality over ar- 
eas of a few /im^ (see Ref. 17). Figure la presents a scan- 
ning electron microscope (SEM) image of a (IIO)YBCO 
film covered with a 3 nm thick (nominal) SRO layer. At 
low thicknesses the SRO layer exhibits island-film topog- 
raphy where the underlying elongated crystallite struc- 
ture of the (IIO)YBCO is clearly apparent in between 
the islands. This island topography enabled us to mea- 
sure, on the same sample, the evolution of the DOS from 
that of a bare (IIO)YBCO surface to the DOS of a fer- 
romagnet in proximity to a (110) surface. Full coverage 
by the SRO layer was obtained at average thicknesses 
of above 10 nm. Below full coverage, the islands fea- 
tured varying thickness even within each island. This is 
apparent in the inset of Fig. lb which focuses on one 
SRO crystallite. It is important to note that in thin 
SRO films of thicknesses lower than 5 unit cells, Tcurie 
decreases substantially^^ Therefore it is expected that 
the exchange field and local magnetization in films that 
feature thickness variations at these ranges will not be 
uniform. Fig. lb depicts a SQUID magnetometer mag- 
netization measurement performed on a 3 nm thick SRO 
layer over-coating a (IIO)YBCO layer. The onset of the 
ferromagnetic transition at ^ 150 K is clearly observed, 
indicating that our SRO layers at these thicknesses are 
ferromagnetic. However, it is also apparent that there is 
a distribution of the value of Tcurie '■ the transition is not 
sharp and takes place over a range of temperatures start- 
ing from 150 K down to 130 K (marked by arrows in Fig. 
lb), indicating that the SRO layer is magnetically inho- 
mogeneous. This inhomogeneity stems from two main 
factors. First is the island topography that yields differ- 
ences between the mean exchange field at the crystallite 
edge and its interior. Second is the thickness variations, 
as discussed above,-!^ The diamagnetism of the YBCO 
sets in at about 88 K. 

The tunneling spectra of the bare (IIO)YBCO films 
featured a pronounced ZBCP with gap-like features all 
over the YBCO crystallites (see Ref. 17 for data on sim- 
ilar films). Interestingly, ABS were also detected in the 
tunneling spectra obtained on the SRO crystallites. Fig- 
ure 2a presents an STM image of a 2 nm thick SRO layer 
overcoating a (IIO)YBCO film where an SRO crystallite 
is apparent. The tunneling spectra presented in Fig. 2b 
were taken at the points marked in the topographic im- 
age. Far enough from the crystallite (point 1), on a bare 
YBCO region, a pronounced centered ZBCP was mea- 
sured (upper black curve), similar to those measured on 
bare (110) films. The middle (red) dl/dV curve was 
obtained at the vicinity of the crystallite edge where the 
thickness of the SRO is lower compared to the center 
of the crystallite (point 2). The amplitude of the peak 
is largely suppressed and it is shifted by ~ 0.8 meV to 
positive energies. The lower (blue) dl/dV curve was ob- 
tained on the SRO crystallite (point 3). Here, the ZBCP 



FIG. 2: (color online) (a) 80 nm^ STM image of a 2 nm t; 
SRO layer on a (IIO)YBCO film showing an SRO cryf 
lite (yellow brighter area), (b) Tunneling spectra (vertic 
shifted for clarity) obtained at the points marked in (a) 



is shifted to +3 mcV and the peak appears to be brc 
ened. In fact, our overall accumulated data (see bel 
implies that this shifted and broadened peak refled 
split ZBCP in which the negative energy peak (in 
case) is largely suppressed and smeared (due to tL' 
and life-time broadening effects). 

The typical tunneling spectra that were obtained 
the SRO crystallites featured a pronounced split of 
ZBCPs with a varying degree of imbalance between 
positive and negative peak heights. This imbala 
ranges from a nearly symmetrical, double-peak struct 
to an asymmetrical split in which one of the peaks is f 
suppressed. Interestingly, the suppression of both 
positive and negative peaks were observed with rou| 
the same occurrence (even on the same sample). M( 
over, the amplitude of the split (energy shift of A. 
also featured pronounced changes ranging between 
than 1 meV and a maximum energy shift of 5 n 
(much larger than the known split amplitude due to 
admixture of a sub-dominant OP at optimal dopin — 
We note that a wide range of both the degree of im 
ance and the energy shift were observed on all sam' 
with average SRO thickness up to 9 nm. Figure 3 sh 
the range of peak imbalance and the suppression of i 
ative as well as positive peaks. Here, the dl/dV cu] 
were all obtained on the same sample (but in diffei 
areas), with an SRO thickness of 3 nm. The peak; 
Fig. 3a (with shifts of up to ~ 3.5 meV) exhibit a smu,^^, 
but pronounced, imbalance between peak heights where 
the positive peak is suppressed relative to the negative 
peak. The imbalance is larger in Fig. 3b in which the 
positive energy peak seems to be fully suppressed and 
the negative peak is found at -3.1 meV. The degree of 
imbalance in Fig. 3c is roughly the same as in 3a but 
here the negative energy peak is suppressed. In Fig. 3d 
the suppression is, again, at the negative side and only 
a trace of the negative peak is still detectable, yet more 
than in 3b. 

The range of the energy shift is demonstrated in Fig. 
4. Here, the dl/dV curves were acquired sequentially at 
equal steps along a 200 nm long line taken on a 2 nm 
thick SRO layer. The ZBCP evolves from a centered 
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FIG. 4: Tunneling spectra obtained sequentially at fixed steps 
along a 200 nm line on a 2 nm thick SRO layer overcoating 
a (IIO)YBCO film demonstrating the range of energy shift. 
The curves are shifted vertically for clarity. 



ZBCP (bottom curve) into a shifted ZBCP where the 
energy shift increases up to a value of 4.5 meV (upper 
curve). This variation of the energy shift occurs on a 
single SRO island (as demonstrated also in Fig. 2). It is 
important to note that there is no clear correlation be- 
tween the magnitude of the split (energy shift) and the 
degree of imbalance between peak heights (e. g., com- 
pare figure 3a and 3b), suggesting that these two effects 
have a different origin. 

We shall first address the mechanism by which ABS 
can penetrate the F layer. The formation of ABS at 
the bare nodal (110) surface in YBCO can be under- 
stood within the framework of a quasiclassical insula- 
tor (I)/N/S quantum well mode&, where an effective N 
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layer of thickness results from the pair breaking na- 
ture of nodal surfaces. Here, the current carrying ABS 
correspond to closed quasiparticle trajectories consisting 
of ARs at the S/N interface and normal reflections at 
the free interface. A metallic N layer, that is in good 
electrical contact with the S, serves as an extension of 
the above-mentioned inherent N region at the nodal sur- 
face and the INS model is still valid. The picture is 
more intricate if the N layer is replaced with a spin po- 
larized F layer. The AR is suppressed due to the im- 
balance between up and down spin populations. More- 
over, the wave- vector mismatch between quasiparticles 
with up and down spins in F leads to breaking of the 
retro-reflectivity of the AR process. — This means that af- 
ter two consecutive Andreev and normal reflections the 
quasiparticle will reach the S/F interface at some dis- 
tance, L, from the originally impinging quasiparticle and 
thus a closed trajectory will not be realized. However, as 
long as L < ^s, ABS can still be formed. The above con- 
straints, and the virtual Andreev reflection process2ii2^, 
will result in a decreased spectral weight of the ABS. Nev- 
ertheless, the penetration of ABS into an F layer is still 
possible and is observed here by the pronounced ZBCP 
we measured in the DOS of the F layer. The reduced 
spectral weight is also demonstrated in the above data 
(figures 2 and 4). The ABS will survive in the F layer as 
long as its thickness is smaller then the mean free path 
and phase coherence is maintained between the electron- 
like and hole-like quasiparticles along their trajectories. 
Our results show a clear penetration of the ABS up to 
thicknesses of 9 nm. The mean free path in our SRO is- 
lands could be significantly smaller than 50 nm, the value 
measured in single-crystals at 1.8 K.— Phase coherence in 
F is maintained on the length-scale of the the F-coherence 
length, which is estimated to be 3 nm in the clean limitj^ 
Therefore, one should not expect penetration to distances 
larger than ^ 9 nm. We emphasize that the measured 
ZBCP was continuous all over the SRO crystallites and 
was not localized along narrow and elongated strips. For 
that reason, the observed penetration of ABS into the 
SRO layer can not be accounted for by effects related 
to magnetic domain walls, as was the case in our previ- 
ous work on (110)YBCO/SRO (antinodal-YBCO/SRO) 
bilayers.^** In that study, the OP penetrated the F layer 
to distances above 26 nm, an order of magnitude larger 
than the F-coherence length. This penetration, however, 
took place only along narrow and elongated strips, sep- 
arated by at least 200 nm, consistent with the known 
magnetic domain wall structure in SrRuOs . This behav- 
ior was attributed to crossed Andreev refiections, taking 
place in the vicinity of the magnetic domain walls 

The observed imbalance between peak heights is a con- 
sequence of the spin polarization inside the F layer: the 
peak that corresponds to the minority spins in F is sup- 
pressed with respect to that corresponding to the ma- 
jority spins. However, the exchange field in F is not ex- 
pected to split the penetrating ABS, as it does not split 
the coherence peaks of the observed^ii^ induced gaps in 



S. We therefore suggest that the splitting of the ABS 
takes place at the interface or inside S. Our observa- 
tion that the negative and positive energy peaks are sup- 
pressed with the same occurrence supports this conjec- 
ture. When the spin degeneracy of the ABS is lifted 
and the ZBCP is split, the peak that corresponds to the 
minority (majority) spins is shifted to positive (nega- 
tive) energies. Consequently, if the splitting occurs in- 
side F, then the spin polarization should always suppress 
the peak that corresponds to the minority spins, i. e., 
the suppression is always expected to occur in the same 
side regardless of the direction of magnetization in F. 
Which of the sides (positive or negative) will be sup- 
pressed in F is determined by the sign of its spin polar- 
ization, which is known to be negative in SRO (i.e., the 
majority spin at the Fermi surface is antiparallel to the 
bulk magnetization).— Therefore, only the negative peak 
is expected to be suppressed in that case, in contrast to 
our observation. This strongly indicates that the split- 
ting cannot occur inside the F islands. Moreover, the lack 
of correlation between the degree of splitting and imbal- 
ance provides further evidence that the splitting does not 
occur inside F. We believe that the observed splitting is 
a consequence of an induced magnetization inside S that 
penetrates it to a distance of ^s- As noted above, this 
induced magnetization can be parallel or antiparallel to 
that in F, depending on local changes in the thickness 
of the F layer and the interface transparency] i i2,i3, 14 
Therefore, both negative and positive sides can be sup- 
pressed (on the same sample), depending on the relative 
direction of magnetization with respect to that of F. 

We now turn to discuss the magnitude of the split. 
The relevant energy scale is the exchange energy, Eex ~ 
KsTcurie, which sets the upper limit for the value of the 
energy shift. In SRO, E^x ~ 13meV and would corre- 
spond to a maximum energy shift of ^ 6.5 meV to each 
side. The induced magnetization should be weaker and 
hence the lower values we typically observe. As demon- 
strated in figures 2 and 4, the observed energy shift can 
change at different lateral locations on the same F layer, 
even on the same crystallite. This can be explained by 
local changes in thickness of the F layer and/or the in- 
terface transparency that cause a non-uniformity in the 
induced magnetization. In general, the energy shift (or 
split amplitude) become larger as we move towards the 
center of the island, where the thickness is usually larger 
compared to the rims. This is consistent with stronger 
magnetization and consequently stronger induced mag- 
netization in S, yielding larger splits (shifts) at the YBCO 
side of the interface. The magnetic non-uniformity of the 
SRO layer, manifested in the distribution of the SRO's 
Tcurie in Fig. lb, supports this assumption. The junc- 
tion transparency, on the other hand, is not expected 
to change within an island over such distances. So, it 
seems that the local magnetization affects our spectra, 
although the interface transparency and its variations 
(mainly from island to island) may also play an impor- 
tant role. 
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A theoretical model that is seemingly related 
to our experiment was introduced by Kashiwaya 
et al., who studied the conductance properties 
of fcrromagnet/ferromagnet-insulator/d-wave supercon- 
ductor junctions. Their calculated I-V characteristics 
exhibit (in the case of the nodal S surface) spectral fea- 
tures similar to some of those observed in our measure- 
ments, namely, an asymmetric split of the ZBCP, where 
the peak at positive energy is suppressed. As in our 
case, the polarization of the F electrode suppresses the 
peak that corresponds to the minority spins. However, in 
their model the splitting of the ZBCP is a consequence of 
'spin filtering', namely, different effective barrier heights 
that are felt by quasi-particles with opposite spins that 
tunnel through the ferromagnetic insulator.— In spite of 
the fact that we could achieve good fits to our spectra 
using this model, we do not believe that it applies to 
our experiment, as explained below. Very low F/S in- 
terface transparencies are needed in order to reproduce 
our spectral features by this model, even lower than the 
typical values corresponding to the tunnel junction be- 
tween the STM tip and the sample {Z ^ 5 within the 
BTK formalism--). Such low transparencies are not con- 
sistent, to say the least, with our previous STS studies 
of SRO/YBCO bilayers^i and magnetoresistance inves- 
tigations in YBCO/SRO/YBCO trilayers^ where prox- 
imity effects were observed, implying highly transparent 
(low Z of less than 1) interfaces (that also allow for in- 
duced magnetization in the YBCO layer) . Consequently, 
in our measurements transport is dominated by the tip- 
sample tunnel junction and therefore our spectra reflect 



the DOS at the SRO surface and thus monitor the pen- 
etration of the ABS into the SRO layer. In contrast, the 
I-V spectra calculated by the Kashiwaya et al. model 
reflect the properties of the S/F interface, in particular 
the effect of spin flltering. Wc also like to point out that 
the spin flltering scenario cannot naturally account, in 
our bilayers, for those spectra in which the negative peak 
is suppressed. For this to occur, one has to impose the 
existence of an insulating ferromagnetic layer (between 
the YBCO and SRO layers) having spin polarization op- 
posite to that of the bulk SRO, which is unreahstic. 



IV. SUMMARY 

In summary, our scanning tunneling spectroscopy 
measurements of SrRuO3/{110)YBa2Cu3O7-s bilayers 
clearly reveal a penetration of the ABS into the F layer 
and provide evidence for the predicted splitting of ABS 
at the S/F interfaces due to an effect of induced magne- 
tization at the S side. 
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